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NONLINEAR EQUATIONS OF MOTION FOR CANTILEVER ROTOR BLADES IN HOVER 
WITH PITCH LINK FLEXIBILITY, TWIST, PRECONE, DROOP, 

SWEEP, TORQUE OFFSET, AND BLADE ROOT OFFSET 
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and 

Ames Directorate 

U.S. Army Air Mobility R&D Laboratory 
SUMMARY 

Nonlinear equations of motion for a cantilever rotor blade are derived 
for the hovering flight condition. The blade is assumed to have twist, pre- 
cone, droop, sweep, torque offset and blade root offset, and the elastic axis 
and the axes of center of mass, tension, and aerodynamic center coincident at 
the quarter chord. The blade is cantilevered in bending, but has a torsi nal 
root spring to simulate pitch-link flexibility. Aerodynamic forces acting on 
the blade are derived from strip theory based on quasi-steady two-dimensional 
airfoil theory. The equations are hybrid, consisting of one integro- 
differentlal equation for root torsion and three integro-partial differential 
equations for flatwise and chordwise bending and elastic torsion. The equa- 
tions are specialized for a uniform blade and reduced to nonlinear ordinary 
differential equations by Galerkin's method. They are linearized for small 
perturbation motions about the equilibrium operating condition. Modal analy- 
sis leads to formulation of a standard eigenvalue problem where the elements 
of the stability matrix depend on the solution of the equilibrium equations. 
Two different forms of the root torsion equation are derived that yield vir- 
tually identical results. This provides a reasonable check for the accuracy 
of the equations. * 

INTRODUCTION 

The general problem of helicopter aeroelastic stability involves coupling 
between the motion of the individual blades and coupling between the rotor and 
the body of the helicopter. The complexity of the general problem poses a con 
slderable challenge to the analyst, both in developing an analytical model of 
the system and in understand lag its physical behavior. An important part of 
the general rotor-body dynamic system is the single blade rotating about an 
axis fixed in space. For many problems of practical interest, blade-to-blade 
and rotor-body couplings are not significant and the analysis of a single 
rotor blade constitutes an important problem by itself. Even when coupling 
with other blades and the body is significant, single blade behavior usually 
remains recognizable and can be helpful in understanding the behavior of the 
more complete system. For this reason, the dynamics of a single blade forms 
an important fundamental building block in the study of helicopter dynamics. 

*Ames Directorate, U.S. Army Air Mobility R&D Laboratory 
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Helicopter rotors with cantilever blades are commonly termed "hinge less 
rotors." In contrast with the more conventional articulated rotor, the canti- 
lever blades of the hingeless rotor are attached directly to the hub without 
flap or lead-lag hinges. This configuration reduces mechanical complexity and 
improves helicopter flying qualities by increasing rotor control power and 
angular rate damping. The lack of hinge articulation also alters the struc- 
tural characteristics of the rotor blade and can signif icantly influence 
aeroelastic stability. 

Aeroelastic instability is possible because of the structural coupling 
between bending and torsion deflections of cantilever blades. This type of 
instability is usually characterized by coupled flap bending, lead-lag bending 
and torsion deflections with a frequency near the lead-lag bending natural fre- 
quency. The structural coupling of cantilever blades is significantly depend- 
ent on the specific configuration parameters of the rotor blade, and the 
magnitude and variability of this coupling make the analysis of cantilever 
rotor blades an important and complex subject. 

A comprehensive study of hingeless rotor stability is a formidable task 
because of the many important configuration parameters. The cantilever blade 
structure treated here is shown in figure 1. The elastic blade can be rotated 
about the pitch change bearing by vertical movement of the pitch link from the 
swashplate controls. Pitch-link flexibility, represented by a spring element, 
will permit, rigid body pitching motion of the blade (i.e., root torsion). Cer- 
tain small offsets of the blade axis are often provided to reduce steady blade- 
bending stresses, to improve rotorcraft flying qualities, or to enhance rotor 
blade aeroelastic stability. Five jf these offsets are considered in this 
report: precone, droop, sweep, torque offset, and blade root offset. Precone 

is the inclination of the pitch change bearing with respect to the plane of 
rotation (positive upward). Droop is an inclination (positive downward at 
zero pitch angle) of the blade segment outboard of the pitch change bearing. 
Sweep is a rotation of the blade in the plane of rotation (at zero pitch angle) 
about the blade root, positive in the direction of blade rotation. Torque 
offset is a lateral shift of the blade in the plane of rotation (positive in 
the direction of blade rotation). The blade root offset is the distance 
between the center of rotation and the root of the blade. These parameters 
are illustrated in figure 2. 

The equations of motion developed in this report are an extension of the 
previous nonlinear equations of Hodges and Dowell (ref. 1) which treated a 
general nonuniform, twisted, torsionally elastic cantilever blade. The configu- 
ration of reference 1 included precone as well as chordwise offsets of the 
blade section mass center, tension center and elastic axes. The present con- 
figuration does not include offsets between the chordwise axes but does 
include the additional configuration parameters: droop, sweep, pretwist, 

torque offset, hub offset, and pitch-link flexibility. For additional back- 
ground discussion of the aeroelastic stability of hingeless rotor blades, and 
a discussion of recent pertinent research, the reader is referred to 
reference 2. 

The equations of motion are derived from Hamilton's principle. First, 
the structural terms are adapted from reference 1; the inertial terms are then 
formulated from the kinetic energy. Blade aerodynamic loads are developed in 
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a way similar to that in reference 2 using strip theory and quasi-steady 
two-dimensional airfoil theory. The resulting hybrid equations are then spe- 
cialized for a uniform blade. The three nonlinear integro-partial differential 
equations and one integro-dif t erential equation are transformed into 3N ♦ 1 
nonlinear ordinary differential equations by Galerkin's method, where N is 
the total number of mode shapes for each of the elastic torsion, lead-lag, and 
flap deflections. When these equations are linearized about the equilibrium 
operating condition, two sets of equations result. The equilibrium deflection 
is specified by 3N + 1 nonlinear algebraic equations and the stability of 
small perturbation motions about equilibrium is determined by a set of 3N + 1 
homogeneous linear ordinary differential equations with constant coefficients 
that depend on the solution of the equilibrium equations. Preliminary calcula- 
tions fo" stability are presented in reference 3 which treats the effects o£ 
twist, ^icch-link flexibility, precone, and droop. 

Technical discussions with Dr. Robert A. Ormiston and checking of the 
entire derivation by Dr. Donald L. Kunz are gratefully acknowledged. 

DERIVATION OF THE EQUATIONS OF MOTION 

In this section we consider a blade for which the structural, inertial, 
and aerodynamic properties may be treated independently. Here we assume that 
the cross section structure 1 and inertial properties are doubly symmetric with 
respect to the blade cross section principal axes. The blade, shown schemati- 
cally in figure 3, rotates at constant angular velocity 0 about the k r axis 
fixed in space. The derivation of the equations of motion is fundamentally the 
same regardless of the assumption of double symmetry of the blade structural 
cross section. The aerodynamic center of the blade is assumed to be at the 
quarter chord and coincident with the blade elastic axis. The deformed blade 
position may be described in terms of several coordinate system rotations, the 
sequence of which is precone 8p C , pitch (0 O prescribed and *(t) restrained 
by the pitch-link stiffness k$, as shown in fig. 4), sweep C 8 » droop Bj , 
and the blade elastic deformations u, v, w, and <). The blade bending deflec- 
tion coordinates v and w are defined parallel to and fixed to tlje blade 
principal axes at the blade root as pictured in figure 5 (along j and k in 
fig. 3). Also shown in figure 5 are the axial deflection u and the blade 
elastic torsion deformation $. This notation is a slight modification from 
that of reference 1 where v and w were parallel to the horizontal and verti- 
cal planes, respectively. The pretent scheme is more convenient when assumed 
modes are used since the vertical and horizontal deflections are functions of 
precone, droop, sweep, torque offset, and the total root pitch angle (0 O + $) 
as well as the blade elastic deflections. 

Structural Terms 

In this section the structural terms are written for a blade with arbi- 
trary radial distribution of stiffness properties. These terms are taken 
directly from the final equations of reference 1. The parameters e^, Bi*, 
and C}*, defined in reference 1, are set equal to zero because of the double 
symmetry of the cross section assumed here. These terms were derived from 
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integrating the variation of the strain energy by parts anti collecting the 
coefficients of Su, 6v, 6w, and 6>. For the coordinate system of this report 
we must let d in the equations of reference 1 represent the bade pretwist 
angle 0 pc , defined to be zero at the blade root. We also assume that 6 pt is 
a small angle. 0(c), where c is of the order of magnitude of the bending 
slopes. The structural terms for the present configuration, neglecting 
warping rigidity Ci, are thus 


<5u term: 

-T* 

6v terms: 

-(Tv’)’ +{[EV - ( Ey -EI v Q(e pt + ») 2 W + (EI Z » -EI y .)(e pt + *)w"}" 
6w terms: 

-(Tw'V +{(EI Z » -SI y .)(0 pt + *)v" + [EI y . + (EI X * -EI y .)(8 pt + *) 2 ]w H }" 
6$ terms: 

-k A 2 [T(6 p t-»- ♦)']*- (GJ*' ) 1 + (EI 2 . - EI y . ) [ (w" 2 - v ,,2 > '0 pt + + v'V] 


( 1 ) 


where T £ EA[u' + (v' /2) + (w' /2)]. The underlined terras in equations (1) 
and below are terms of O(e^) associated with $ and 6$. These terms are 

small and are retained only for the sake of completeness. The only structural 

term in the root torsion equation may be derived from the potential energy of 
the spring (l/2)k$* 2 . Thus, 

6# term: (2) 

The boundary conditions, strictly from structural considerations are 

x • 0: v-w-v' -w' ■ $ ■ 0 l 

I (3) 

x ■ tl v" ■ w" ■ v’" * w'" ■ "0 1 

Another boundary condition at the blade root is that of torsion moment 
equilibrium. This requires knowledge of the various blade axes systems . The 
rotations of axes are described fully by the following transformations: 
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where C and S are the cosine and sine, respectively, of the subscripted 
angles; [T] is the deformed blade coordinate transformation of reference 1 
with 0 - 0p t ; and for small angles and deflections the T vectors are 
approximately along the blade, the j vectors are approximately horizontal, 
the k vectors are approximately vertical; T * j • k. For determining the 
boundary conditions at the blade root, we note that the pitch-bearing axis is 
along T p . Thus, M ■ M^ip. At the root of the blade, however, 

+ MyJ + Myi< and this, of course, must be balanced by the spring 
moment M$Tp. Here M|, My, and My are twisting and bending moments due to 
w, and v deflections, respectively. Thus, 

M, - M,I • t p + • T p + M,,it • 1,, 

' Vi. C 8 d ' MvSc. + V^S,,, 

* ", - C.k, + e d K, (5) 

Therefore, in terms of deformation quantities 

k** * GJ^ , ^+T(O)k A 2 l0p t (O)+^ , (O)]+C a RI y *w"(O) + e d EI 2 .v"(O) (6) 

This is the torsion moment equilibrium boundary condition assuming that 
0pt(O) and 0p t is a small angle. Equation (6) may be used as a replace- 
ment for the root torsion equation to be developed below. Numerical results 
are virtually identical regardless of which of the equations are used. All 
structural quantities in the above equations are defined the same as those of 
reference 1. 
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•'iAVh'^V^ ■* 


x(e pt + *)’' 

n V;, •• * • 


where x Is the axial coordinate measured from the root of the undefo 
beam to the tip; u, v, w, and 4 are elastic displacements; X is the warp 
function (X ■ 0 at the elastic axis); and n end £ are the beam cross- 


section principal axes centered at the elastic axis. The variation of the 
tem kinetic energy Is then 


- / 1 ff p$ • 6$ dr uc 


We now write the velocity in the l r , jf r , ltj. system 

5 - [-o. 0 - (n + Se^ijA,, + c epc Ai - s Spc Ai]t r 


PC 7 


+ ^ c 0r, c ^ + + s B oc A xl l S: 




' 


r 5* 


m 






Inertial T 

In this section the inertial terms are developed fr- 

trary radial distribution of mass properties. For the 

tural terms we were able to utilize chose derived previously 
This is not feasible wich the luertial terms because of the i 
eteri and the new coordinate system. Thus, we will derive the 
independently from the kinetic energy following the method of i 

The first step is to write down expressions for the 
point In the blade. Symbolically, this la quite simple 


^ • n ^r * <*i^r ♦ «pJr ♦ *it + XiJ ♦ *i^> 
+$t p * (*i 1 + yjT + *i^) + *► yjl ♦ *i 


- v 

■ ' W tv / 


where xj v yj, and z l are the coordinates of a generic point in the 

along the T, J, and lc vectors measured from the root of the beam 
pitch chang. :» cring). Here e 0 16 the torque offset, e; Is the 
set, a J raining unit vectors are given In equations (4). The f 

of equation u) ia the contribution from blade rotation. The flgft tei 
second line is the contribution from root torsion. The final thre#tt 
the contributions from blade deformation rates. Expressions for xj , 
Zi may be obtained from the deformed blade coordinate transformation 
(ref. 1). 




where 


A x * x i c c s c e d - yi s ; 8 + z i c c 8 s 6 d 

Ay - x l( c e 0 +* s t; g C8 d + s 6 0 +^ 6d ) + - Sq^C^) 

A z " X 1 ( S 9 0 +« S i; 8 C B d “ C «^ S B d ) + yi S 0 o +® c C s + 7 ‘1 ( S 9 0 +^ S C 8 S B d + c e 0 +* c S d ) 

A X 3 X l C C n C 8 d " yi 3 Cg + i l C Cg S B d 

^ " Xl ^ C -o + * S ^s C B d + S 9o+* S B d ) *yi c e 0 +* c c 8 + i l( c e 0 +* S Cg S B d ■ s 0 o + * c B d ) 

A £ 3 x i( s e o +*s Cg c 6d -ce o+ ^se d ) + y 1 Se o+ *c C8 + 2 1 (s 9o+(> s C8 s Bd ^C0 o+ *c 6d ) 

Note that with this <5 tort hand notation L. - Aa - *A Z ; «A- - A*„ + A v 6*, etc. 
Thus, the variation of the velocity is simply y y 

«V - [-(n + Sep^XAfiy-Ajfit)^^^ 

■ S 8p c ( A 6 z + A ^ 6 * + Ay«*)] i r 

+ [ fiC B pc A 5x " (ns 6pc + i)(A 6z + Ay««) - A z 6* + A^ - A*«*] } r 

+ ( C 6pc [V* + ( A «y “ A z <5*)$ + A si + A^4>] + S 6pc A 5 *|^ r (12) 

Substitution of equations (10) and (12) into equation (9) yields (with a 
change of sign and integration by parts to accommodate Hamilton’s Principle) 

™ ■ So //a c - n2 =l<V - n^SpeA* - 2!JC apc (A. - iAp) + aJa 4x 

+ (j-n 2 (e 0 + Ay) + 2TiCg pc A x - 2USg pc (A| + $Ay) - (#A Z + $Aj + * 2 Ay) + AyJ a 4 y 

+ [ fi2s 8p C c B pc A x + « 2 =lS 8pc - n2S B pc A* + 2aS 8 pc (Aj - «A Z ) 

+(tA y + tAy - t 2 A z ) tAj] A 4z +[n 2 e 0 A z + a 2 S 6pc C 8pc A x Ay + n2. 1 S Bp<; Ay 
4 fi2c 0pc A y A z + 2J % C (AyAy + A z Af ) - 2nc 0 pc A z A£ + li (A y A^ + A z Aj) 

+AyA^ - A z Ay + $(A y 2 +A z 2 )Jfi*Jdn d? dx ( 13 ) 
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Equation (13) contains an exact representation of the inertial terms. It is 
at this point that one must make some assumptions regarding an ordering scheme 
in order to have a tractable set of equations. The same ordering scheme as 
that of reference l is assumed here, with the addition of Bj, C s , e 0 /S, ej/", 
and ♦ all being 0(e). Here e is a small parameter of the order of magni- 
tude of the bending slopes. We assume that 6.f is truncated so that terms 
0(e 2 ) are neglected with respect to unity. This scheme is not perfect ond, 
hence, is not followed rigidly. It does, however, greatly simplify equa- 
tion (13) and yields a set of nonlinear equations with only the most important 
nonlinear terms Included. 


Even though they are very small terms it has been suggested (ref. 4) 
that higher-order terms should be Included in the root torsion equation. We 
will follow that suggestion for the sake of completeness. Unlike reference 4, 
however, corresponding terms in the other equations will also be retained to 
insure that modal mass and stiffness matrices are symmetric and that the modal 
gyroscopic matrix is antisymmetric (in vacuo ) . The analogous small terms 
associated with the blade elastic torsion equation, underlined in equation (1), 
are also retained on this basis - strictly for the benefit of the doubt in 
both cases. 

We first consider terras in equation (13) which multiply 6xj, dy^, and 
6zi, consecutively. Since 6xj, 6y lt and contain du, dv, dw, and 6$, 

the Inertial terms from these respective partial differential equations may be 
recovered. We need all terras up through 0(c) for 6xj and up through 0(e 2 ) 
for terms involving ♦. Thus, 

6xj terms: 

J Q l JS k -p[fi 2(A x + ei) + 2fi(A£ - *A z )]dn de dx (14) 

We need all terms up through 0(e 2 ) for dyj and 6zi, and up through 0(e 3 ) 
for terms involving $. Thus, 

6y 1 terms: 

CIA pKc B (Ax +e l> + 2 i! « 8 < A ?- iA z> 


+ [-0 2 (e o + A y ) - 2n8p C (A £ + $A y ) + 2fiA* 



-$A Z + Ay]Ce o+ * + [^BpctAx + eO + 2fiB pc (Ay - *A Z ) 


+#Ay + A z ]Se 0 +$}dn dc dx 


(15) 


&z\ terms: 


n pj-fl^dtAx + ei) - 2fiB d (A^- *A Z ) 

+ [n 2 (e 0 + Ay) +2J28 pc (A i + *A y ) - 2«A* + $A Z - AyJSg^ 

+ [n 2 B pc (A x + ei ) + 2nB pc (A^- iAj) 

+^Ay + A z ]CQ o +(j»|dn dC dx (16) 


The 3$ contribution as it appears in equation (13) is unchanged for this 
step in the analysis. Rate product terms such as $ 2 Ay or $A£ are neglected 
since they do not contribute in a linearized stability analysis. 

Secoud, we may substitute equations (11) into equations (14)-(16) and in 
the part of equation (13) to produce 

6xj terms: 

~ P [A 2 (x + ex ) + 2ft(y 1 C0 o + ( j> - zi s 0 o +#) 

-2fii(y 1 S0 o + z 1 C 0o ) - 2fixi(i: 8 S 9o - 6 d C 0o )]dn dc dx (17) 

6yj terms: 

f 0 l Jf A P[« 2 ei(B pc Se o+ ^+5 s ) -n 2 e o C 0o+$ + fi 2 x(B pc + C s S eo+$ - B d C0 o+ *)Se o+ $ 
+2flxiC 0o +(t - 2Qi\ (Bp C + C s S 0o +$ - 8 d C 0o +$) - 2f2xC g B^ < t 
-2GB pc y|*- 2fiC s *(y 1 S 0o + z 1 C 0o ) - Q 2 y 1 Ce o +« 

+n 2 z 1 S 0o+4 C0 o+ ^ - (zx - B d x)* + y 1 ]dn dc dx (18) 
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6z i terms: 

ai p [ fi 2 e! (6p C Ce o+$ - 8 d ) + n 2 e 0 SQ o+ * + n 2 x(8 pc + c 8 se o +# - B d Ce o +$)Ce 0 +$ 

-2nxiS0 o+ ^ + 2fty j (8 pc + C 8 S 0 O +^ “ 8d c 0 o +*) + 20 xB d Bj* 

-2nB nc z 1 i + 2QB d *(y 1 s 6o >z 1 c 0o ) + fi 2 y 1 Se o +$Ce o+ * 

-ft 2 z 1 Se o + 4 »+ (y J +C s x)i + z 1 ]dn d£ dx (19) 

6$ terms: 

/o7/a p [^ 2 xe o (C a Se 0+ * + B d Ce o+< j) + J2 2 e o (y 1 S0 o+4 + EiCe o +$) 

+n 2 x“" ^B pc + C a S0 o +$ - B d C0 o +<p) (C 8 C0 o+( j + B d S0 o+ $) + fl 2 xeiBpcCj 
+fi 2 x(Bp C + SsS0 o +$ ~ B d C0 o +^) (yiC0 o +$ - ZiS0 o +$) + Q 2 ejBp C (yjC0 o - ziS 0 q ) 
+ft 2 x(c s C0 o +$ + B d S0 o 4.^)(y 1 S Qo+<t + z 1 C0 o+( f) + 0 2 (y 1 2 - zi 2 ) S 9o +<j>C0 o +$ 
+n 2 yiZi c 2(0 o+4 ) + 2 nB pc (yiyi + zizi) + 2 «xB 1 (^ 8 y 1 - 8 d *i) 

-2fix(C s S 0o - B d C 0o )x 1 - 2fi( yi S 0o + z 1 C 6o )x 1 

+2«(y 1 S 0o + z\C 6q ) (c s yj - B d Zi) + ( yi + Cgx)^ - (z\ - 6 d x)yj 

+x 2 $(c 8 2 + B d 2 ) + 2x*(c s y 1 - B d Z!> + (y x 2 + Z! 2 )*]dn d£ dx (20) 


where 8 X - B pc + C s S 9o - B d C 0o ; - C 8 C 0o + B d S 0<) . Inequations (17)-(20), the 
expressions have been truncated to the appropriate order of magnitude as pre- 
scribed by the ordering scheme discussed above. Finally we must write 
explicit expressions for xi , yi, and z\ and use these expressions to formu- 
late the Inertial terms explicitly in terms of u, v, w, <fr, and ♦. From 
reference 1 and equation (8) 
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X} ■ x + u - X(0 pt + $) ' - n[v' cos(0p t + ♦> + w' 8in(6p t + $) ] 

-C[-v' 8in(0 pt + ♦) +w' cos(0 pt + ♦)] 
yj « v + n cos(0 pt + $)-C sin(0 pt +<|>) 

zi - w + n 9in(0p t + <t) + C coe(0 pt + <J>) ( 21 ^ 

The 6u terms may be formed directly from the iSxj terms 
<$u terms: -mft 2 (x + e! ) + 2mfixi(C s S 0o - B d C 0o ) 

-2mfl(vC 0o+<t -wS 8o+ *) +2mn(vS 9o + wC 9o )* (22) 

where m - // A p dn d C and ff A pn dn d; - ff A PZ do dC - 0 for coincident 

elastic and mass center axes. The 6v and 6w terms may arise from both x l» 

5y lt and Sz±. The terms from <$xj are small, however, and are underlined 

below 
5v terms: 

n^ 2 ei (6 pc Se 0 +« + C 8 ) - mft 2 e o C 0o +* + raft 2 x(e pc + ; s S 6o +$ - Bd c 0 o +*> s 0 o +$ 
-mJ3 2 vCe o+$ + mfl 2 wS 9o+ ^C 6o+ ^ + 2mftuC 0o+ $ - 2mflwC8 pc + S s S 0o+ $ - Bd c 0 o +$) 

-2milxC s 8 j ♦ - 2mft0 pc vi - 2mfiC s ( vS 0q + wC 0q ) ♦ 

-m(w - B(jx) $ + mv + 2flS 0o (mk^ 2 $) ' + 2flis 0o (mk m2 ) ' (23) 

6w terms: 

mft 2 ei (8p C c e 0 +$ " ®d^ +n^ 2 eo s 0 o +* + mfl2x ( 6 p C + ? 8 S0 o^“ 6 d c 0 o +« )C e o +* 
+mfl 2 vS 0o +^C0 o +^ - mfl 2 wS0 o +^ - 2mfluS 0o+ $ + 2mQv(0p C + C 8 s 0 o +$ - 0 d C 0 o +*) 
+2mflxB d B 1 * - 2mfiB pc w*+ 2mfl8 d (vS 0o + wC 0o )* 

4m(v + ? 8 x)* + mw+ 2JlC 0o (rak4 1 ^) ' + 2n*C 0o (mlc^ 1 ) 1 (24) 

Since the 6v and 5w terms from 5xj arise from integration of 6& by parts, 
there are additional tenu3 evaluated at the boundary that do not necessarily 
satisfy equations (3). These terms are 

-2mft(* + *)O4 2 S 0o «v + k^Ce^w)!* (25) 
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These terms must be included in the analysis when Galerkin's method is applied 
in the solution of the equations of motion. 

The 6$ terms arise from both 6y j and 6z^ since contributions from 
6xj are negligibly small. The contributing terms from 6y 1 and 6zj are 

y.fiy! + z\&z l + (2flxi - n ? -yiC 0o+$ + n 2 ZiS 0o +$) (6y 1 C 0o . f(t - 6z 1 S 0c+(t ) + (y^zj - z^yj)* 

(26) 


The 6$ terms follow immediately: 

t 

6* terms: mSl 2 (k£ 2 - k^ ) [S 0o C0 o + (6 pt + * + $)C 20o ] 

+2mfi(v , S 0o l4 2 

+mk^($ + ♦) (27) 


Here as above the underlined terms are 0(e 3 ). From equation (20) the 6$ 
terms may be written 

64> terms: 

f o !i \ mfl 2 xe o (i: s S0 r r $- 6 d C0 o+ $) + mfi 2 e o (vS0 o+ $ + wC Qo+ <t) 

-Hnfi 2 x 2 (Bp c + Cs^0 o +$ “ ^d^Qo+l*^ ( ( »s^0 o +4 ®d^0 o ^^ m ^^ e l®pc^l 

+m 0 2 x(B pc + ^ g S 0o + 4 > - B d C0 o+ 4,) (vC 0oT<t - wS 0o+4) ) + mfi 2 e 1 B pc (vC0 o - wS 0 q ) 
+ n'fi 2 x(C s C0 o +<j ) + B d S0 o +,j,) ( vS 0 o +4> + wCq q+ q) + mft 2 (v 2 - w 2 ) S 0 o+( j ) C0 o 4. ) j 
+ rafl2vwC 2( 0o+ ^) +mfi 2 (k 2 2 -kjj) [Se o C0 o + (0 pt + ^ + *)C 2 0 o ] 

+2mft8 pc ( vv + ww) + 2mflx8 j ( C 8 v - B d w) 

-2mOx(C 8 S0 o - B d C 0o )u - 2mn(vS 0o + wC 9q )u 
+ 2mft(v'S0 o km 2 +w'C 0o k^) + 2mfl(vS 0o + wC 0£j ) U 8 v - 8 d w) 

+m(v + C 8 x)w- m(w- B d x)v + mk m (^ + $) 

-Hnx 2 (? s 2 + 8 d 2 )# + 2mx(C s v - B d w)4 + m(v 2 +w 2 )#|dx 


( 28 ) 


The 64> terms are under an Integral because 6$ is not a function of x. 
(Hamilton's principle results in a partial differential equation only when the 
6-quantity is a function of x. ) 


Aerodynamic Terms 

The aerodynamic lift and pitching moment acting on the blade in hover are 
based on Greenberg's extension of Theodorsen's theory (ref. 5) for a two- 
dimensional airfoil undergoing sinusoidal motion in pulsating incompressible 
flow. The rotor blade aerodynamic forces are formulated from strip theory 
and only the velocity component perpendicular to the blade spanwise axis 
(the x'-axis in the deformed blade coordinate system x' , y', z' in fig. 6) 
influences the aerodynamic forces. A quasi-steady approximation of the 
unsteady theory for low reduced frequency k is employed in which the 
Theodorsen function C(k) is taken to be unity. The steady induced inflow for 
the rotor is calculated from classical blade-element momentum theory. These 
simplifying assumptions are judged to be adequate for low frequency (mainly 
determined by the blade bending frequencies) stability analyses of a hovering 
rotor. 

In Greenberg's theory ^ref. 5), a two-dimensional airfoil is assumed to 
be pivoted about an axis which may be distinct, in general, from the aerody- 
namic center axis. The airfoil is pitched at an angle e(t) to the free 
stream flowing at pulsating free-stream velocity V(t). The airfoil is verti- 
cally displaced with velocity h(t) positive downward as shown in figure 6. 

The relations for lift and pitching moment per unit length may be expressed in 
terms of the circulatory and noncirculatory components 


L " Lc + L nc 
M - Me + Mjjq 


(29) 


With the airfoil pivot axis (.analogous to the rotor blade elastic axis) at the 
airfoil quarter chord (the airfoil aerodynamic center) these components are 


l NC “ 
L C - 

Mnc " 
Me ’ 


PepSC 

2 

P^ac 

2 


f (h + Ve + Ve + f) 

v(h + Ve + -y) 



(30) 


/ 
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The Theodorsen function C(k) has been set equal to unity in the circulatory 
lift. It should be noted that e is the angular position of the airfoil with 
respect to space; e and e are the angular velocity and angular acceleration 
of the airfoil. The Instantaneous angle of attack of the airfoil 
a - tan" 1 (Up/U-j) is the angle between the airfoil chord line and the resultant 
fluid velocity U of the airfoil. The airfoil velocity components in the y', 
z* principal axis system are Uj and Up shown in figure 6. It is desirable 
to express the aerodynamic forces and moments in terms of Up and U T . 

Assuming that the angles e and a are small yields 

Up - -h - Ve 
U - /u T 2 + Up 2 - V 

Substitution of h and V from equations (31) into equations (30) yields 


(31) 


l NC 

L C 


p« ac 
■ ' ... ■ 

2 

p_ac 
m USS 

2 



(32) 


Next we consider the total aerodynamic forces in directions parallel and 
perpendicular to the airfoil chord line. The noncirculatory lift is taken to 
act normal to the chord line, and the circulatory lift is taken to act normal 
to the resultant blade velocity U. An aerodynamic profile drag force per 
unit length, based on a constant profile drag coefficient c<j and acting 
parallel to the resultant blade velocity, is Included. 


D 


o m ac do , ^ <5 

■in ~r (u t + °p 2 > 


(33) 


The force components and directions are shown in figure 7. The force compo 
nents T, normal to the airfoil chord line, and S, parallel to the airfoil 
chord line, are therefore 

T ■ Lq cos a + Lfl C + D sin a 
S ■ -Lq sin a - D cos a 


(34) 


From figure 5, 


cos a 
sin a 


UT ' 
U 

Up 

U , 


( 35 ) 


Substitution of equations (32), (33), and (35) into equations (34), with 
c<j 0 /a neglected with respect to unity, yields 
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T - HT [- U P U T + f M - f Up + (f ) 2 E] 

S . ( Up2 . | M . !& u/) 

The expressions for aerodynamic pitching moment components may be written from 
equations (30) and (31) as 

«c “ - ~r (if V 


where U has been approximated by U-j in M(;. The total pitching moment is 
then given by 



The aerodynamic force and moment acting on the blaoe at a point on the 
deformed beam elastic axis (coincident with the blade airfoil section aero- 
dynamic center) are 


dynamic center) are 

t « sj' + Tk' (39) 

M - Mi' (40) 

The virtual displacement and virtual rotation are given, respectively, by 

6q - <5ut + 5v^ + 5wk + 6*i p x [(x + u)t + v^ + wit] (41) 

Sul ■ (42) 

Thus, the total virtual work of the aerodynamic loads is 

\ 

<5W - f * ($ • 6q + M • 6w)dx (43) 

*o 


In order to write the aerodynamic terms in each equation we must express equa- 
tion (43) in terms of u, v, w, <J>, and *. This entails writing Up, U T , and e 
in terms of u, v, w, $, and ♦ . The blade airfoil velocity and rotation are 
simply 

V - -fle 0 l r + ftei + v ± tc r + ut + v} + w£ + (flk- + i>i p ) * [ (x + u)t + vl + wk] \ 

} (44 > 

w ■ flk r + iip + 4 , i’ I 
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where Is the induced Inflow velocity equivalent to that of reference 2 

except that root torsion is Included and the pitch angle is not necessarily a 
small angle. From figure 6 

Up - V • k’ \ 

U T - V • J’ \ (45) 

e • u> • t* / 

Without writing all the details, substitution of equations (36), (38)-(42), 
(44) and (45) into equation (43) yieldj an expression for 6 W of the form 

<$W * (L u 5u + Lyfiv + lv<Sw + M^ 6 $)dx + 6 $ dx (46' 

where L u , Ly, and Ly, are aerodynamic forces per unit length and M^, and M$ 
are aerodynamic moments per unit length. All these quantities appear as 
forcing functions on the right hand side of the equations of motion. We 
assume 0 pt • 0(e) and 0 O » Ofe 1 / 2 ) so that the expressions are not overly 
complicated by small terms. 

L u 2 n 2 x " s e o c 0 o w (A7 > 

Ly - ^ (n 2 (x 2 + 2eix)Se o + n 2 x 2 (e pt + *3S0 o C0 o 

+Vi 2c J o - nxv 1 [s 2 e o + (0 pt + <j>)c 20o ] - ^ n 2 x 2 c* o 

+jo 2 x 2 s 28o + [fi 2 x 2 ( 0 p t + $ + *)- 2 nxv 1 ]c 26o }« 

+ {( 0 pt + *)[«x(S 0 o + *C 0o )-v 1 Ce o ]-2 ^ «xc 0o Jv 

+ {2y, i CQ o - fix[2S 0o + (2* + 0 pt + $)Cq o ]}w> (48) 

Lw * ^ 2 (x 2 + 2e 1 x)C 9 o [S 0 o +(0p t + 0)Ce o ] +fi 2 x 2 C 0Q f* v'w" dx 

-n(x+e 1 )v 1 C 20o - n 2 e 0 x( 6 pc C 2 e o - 8 d C 0o ) + fl 2 f (Bj + w'C 9 o )C 6o 
-fl 2 e 0 xw'Cg o - n 2 xvw'C 0o - n 2 xS 1 vC 0o - fl 2 x 2 8 pc 4 a C 0o + n 2 (x 2 + 2 eix)$C 2 0 o 
+ l n x[Se o + (* + 20 pt + 2*)C 0o ] -v 1 C 0 o }v-fl(x + e 1 )wCe o -|- w 
+ ^ nxC 0o i + [-^- (v+C 8 x)]nxC 0 o *> (49) 
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(50) 


- - nr (»)* *n 2 xc 9o (e, t»'c 8o ) -«l 

V «« <«* - C HT (*(i) 2 t2«*C8o < » + i > + “ 2x < B l +w ' C 9o >1 

+ (f)[f ■ < v+ c jX >] w+ r.-(x J + 2e,x)So 0 (-S 9 o (w - 6 d x) + C 6() <v + C s x) 1 

+fl 2 (x 2 + 2 e 1 x)C 0 o ( 0 pt + ♦) (v+ C s x) - O 2 x 2 ( 0 p t ♦ ♦)S 9 Q Ce o (w - B d x) 
i Cd \ 

-fvi 2 -n 2 x 2 j C 0 Q (w- 6 d x) +n 2 x 2 Ce o (v + ^ g x)jT X v'v" dx 


-ft(x+ei)v l C 20o (v + c s x) +Oxv 1 [S 2 0 o + ( 0 p t + 4 >)C 20o ](w- 0 d x) 
-fl 2 e o x(8 pc C 2 0 o - 8 d C eo )(v-*-C 8 x) ♦ y- ( 0 ! +w , C 0o )Ce o (v + c 9 x) 
-fl 2 xC 0o (e o w’ + vw'C0 o + BjV + x6 pc ; s ) (v + C 8 x) + 2 fl 2 e 1 x*C 26o (v + r JgX ) 
+fl 2 x 2 «>[(v + c 8 x)C 2 0 o - (w- 8 d x)S 2 0 o J - [« 2 x 2 (0 pt + <|> + *) 
- 2 flxv 1 ]^C 20o (w- B d x) + v ({2 OxC 0o - ( 0 pt + O)[nx(S 6o + *C 6o ) 

" v i c 0 o l) (w " s d x) + l nx f s 0 o ♦ <♦ + 29 pt + 2 ^> C 0 O ] “ v i c B 0 }< v + ^s y) ) 
+w<{0x [2Se o + (2* + 0 pt + ♦)C 6o ] • 2v i C 0o [ (w - B d x) 


-0(x+ ei)C 9o (v + (; 8 x)^ + ~ flx(v + 5 g x)C 0o $ 
+ £-^r- (v+C g x)J (v + ; 8 x)OxC0 Q i^dx 


where 


n * ±nR f { / 1+ ^ | s 8„ + M!) + *o(z) + *] c e„| * 1 } 

sgn v t - sgn |s 0£) + [0 pt (|)+* 0 (f) + *] C e 0 j 


be 

irR 


(51) 


(52) 


Equations (1), (2), (22)-(24), (27), (28), and (47) — (51) wher combined as 
structural terms + Inertial terms - aerodynamic terms ■ 0 
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yield the hybrid equations of motion. The boundary conditions are found In 
equations (3) and (6). 


SOI, IT ION OF THE EQUATIONS OF MOTION 


Simplification and Nondlmenalonalizatlon 

The equations are solved by Galerkln's method using mode shapes of a 
nonrotating uniform cantilever beam. For convenience we treat only the case 
with uniform radial distributions of mass and stiffness properties. In gen- 
eral, a flexible beam segment Inboard of the pitch-bearing Is present and non- 
uniformities in bending stiffnesses influence the flap-lag structural coupling. 
In this paper, the inboard beam segment is eliminated entirely by placing the 
pitch change bearing In the hub Itself, thus leaving only the single outboard 
blade segment. However, the effect of the inboard beam segment on the struc- 
tural coupling between flap and lead-lag bending is represented in an approxi- 
mate fashion. Flap-lag structural coupling uopends on the relative stiffness 
of the blade segments inboard and outboard of the pitch-bearing because the 
principal elastic axes of the outboard blade segment rotate through the angle 
3 0 as the blade pitch angle varies, while the inboard segment principal axes 
do ..ot. The resultant effective orientation of principal axes depends on the 
bladi geometry and distribution of bending stiffnesses inboard and outboard of 
the pitch-bearing. Although the variations in the structural coupling signifi- 
cantly influence stability, they are not present in a simple single-segment 
uniform beam treated here. They are difficult to include exactly without 
resorting to a more general blade configuration and a more sophisticated anal- 
ysis. However, an approximate representation of these effects may be intro- 
duced with no increase in complexity. This is accomplished by arbitrarily 
assuming that the averse inclination of the principal elastic axes of a non- 
uniform blade is equal to some fraction of the Inclination of the principal 
axes of a uniform single- segment blade. This entails having the structural 
principal axes inclined at <R0 O rather than 8 0 while the mass and inertial 
* rm8 are unchanged. The factor tf? is called the structural coupling param- 
eter. When • 1, the original equations are retained, but as tf? is reduced 
to zero, the flap-lag structural coupling terms diminirh and eventually vanish. 
Although this is only an approximation of the true effect of flap-lag struc- 
tural coupling, it greatly simplifies the numerical model and does represent 
the type of behavior that would be exhibited by a general nonunlform blade. 

In order for the structural axes to be inclined at <?0 O to the plane of rota- 
tion, however, we must substitute 6pt - (1 - <R)9 0 for 9 in the structural 
terms of reference 1 because of the change in the coordinate system. 

The equations are further simplified by eliminating the 6u equation and 
all terms containing u. This is accomplished by solving the 6u equation 
(consisting of equations (1), (22), and (47) combined as discussed above) for 
T. The resulting expression for T may then be substituted into the other 
equations. An expression for u may be easily obtained using the definition 
of T. Proceeding in this manner we first write the 5u equation 
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-T 1 - mfl 2 (x ♦ e j ) ♦ 2mftx6 ( C s S 9q - e d C 0<) ) 

-2rafl(vCe o+t -wS 0o ^)+2nf2i(vS 9o + wC0 o )+^ Q^Se^w 1 -0 (53) 

For a uniform blade we may nondlmena tonal ize aa follows 

-T • - (x + 5 , > + 2x( C.S, o - 6 d c„ o > ♦ + 2 (5S flo + wc 6o ) ♦ 

-2(5C 9o+ *-4s 9o+ »)+J Se^g^w 1 -0 (54) 


where 


T - 


tnfl 2 ! 2 


4 ’ 


u V 

u " l ’ v “ t ’ 


w 


e l . ip^aei _ lock number 

61 " * ’ Y " ® = (R 4 - ej 4 )/i 4 


Now, integration yields 

1- x 2 


♦ (1 - ^)ei - (1 - x 2 ) (C 8 S 0 - 8 d C e )♦ 


-2 if- 1 (vSe o + wCe o )dx + 2 f_ l ($C 0O ^ * s 9 0 +6> d * 

• 6 S *o C 9o fi * 2 *' d * 


(55) 


Equation (55) for T may now be substituted Into the bending and torsion equa- 
tions retaining terms to the appropriate order of magnitude. In the bending 
equations the contribution of the last term, an aerodynamic term proportional 
to y, is negligibly small. In the torsion equation only the first term is 
retained. 

In order to eliminate u from the equations, we now consider the defini- 
tion of T, equation (If), nondimensionalized with respect to mft 2 & 2 


EA v' 2 w’ 2 \ 

T -3*1* V 1 ~ 


(56) 


Differentiating equation (56) with respect to 6 yields 


u ■ ■ t - v*v* - w’w’ 



(57) 


19 


Typical hlngelesa rotor configurations hava sufficiently large radial stiff- 
ness EA so that nfl 2 t 2 /BA - 0(e 2 ). Thus the first tens of aquation (56) is 
quite a isall whan substituted into the bending and root toraion equations con- 
taining 2mftu and for our purposes it is negligibly small. Therefore* u is 
obtained by integrating equation (57) 


u • “f Q X (v’v' + w'w')dx 


Equation (58) may now be substituted for 6 in the bending and root toraion 
equations. The nondimens lonal equations for v, w, and ♦ now nay be 
written for a uniform beam 


6v equation! 


A 2 v- - (A 2 - Aj)|v n,, s| o - S 25o [(0 pt ♦ ♦)v*J M -»- C 2 ^ o [(9pt ♦♦) 2 v h ] m | 

I S 25q I 

-5"" C 2 J o [(0 pt + ♦)«"]" + S 29(j [(8 pt + ♦) V]" | 

I 

-(tv' ) ' + «i (8p C 8e o +q + C 0 ) - e 0 Ce o ^g + x(0p C ♦ t. a S Q +$ - 8<|Ce +$)Se +$ 

-vCe 0 ^ + wS9 0 ^C eo+# -2Ce o ^ f* (v 1 *’ ♦ w'w’H*- 2w(Bp C + c,s 0 +*-a d c e +*) 


-25c,6j* - 28 pc »* - 2t,(vS, 0 + 5Ce 0 )i - (X- 8 d x)K*5+ 2kJ 2 s 9o i’ 
*%( <**i c 28 0 - ***8 O «0 O )» - {* J S28 0 + I* 2 < 6 pt ♦♦♦♦>- 2 *vi]C 29o |* 


+ 1 2 -r - (e Pt + ♦) »<*„„ + *C eo ) - »iC 9o ]J» 

+ j»[2S So + (28+ 8 pt + pjCjo] - 2 »iC 9 

[ pa 

( 5 2 + 2 *,x)S 9o + 5 2 0 p t S 9 o C 9o + Vc 9 „ - xVi(S 29o + e pt C 29o ) - -f ; ! C , 2 I ( 59 ) 




Wrt.K »r 
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6w equation: 


I S 20 

-v"” ~2 £ -*-c 2 e 0 [(e pt + *)*"]" + s 2 e 0 [(e pt + t) 2 *" 1" 

+A 1 w" M + (A 2 - A ! ) {w :), ’s| o - S 2 $ q [ (0 pt + ♦ )*"]" + c 2$ 0 [( 9 pt + l 

-(tw')' +e 1 (6 pc C 0o+ ^- B d ) +x(8 pc + ; # S0 o+ ^- fJ d C 0o+ ^)C 0o+ ^ + e o Se o +® 

"^ s O o +^C 0o+ ^ - w 80 o +^ + 2S 0o+ * f Q * (*•*' +w , w’)dx + 2v(0 pc +; g S 0o+ ^- B d C 0() ^) 
+2x6 d 8 j ♦ - 20 pc wi + 28 d (vS 0q + wC 0() ) • + (v + C s *) ♦ + & ♦ 2k^ 1 C 0o * ' 

+ ^ ( x 2 + 2ei5 )c e o < t'-x 2c 0 o f* v’w" dx-^ ^ 0o + e o 5w , C 0o + xw , C 0o 

■*2 • 
•»-5e 1 vC 0o - (x 2 + 2e 1 x)$C 20o - (|) (♦ + *) + {v ± C 0o - x[S 0q + (* + 20 pt + 2f)C 0o ]f v 

+ (x + ex)C 0o w + J w-~ xC 0o $- [^- (vtC 8 x)]xC 0o t) 

■ 6 [ (52 + 25lX)c e o (s 9o + 0 pt c e o ) 

-(x + e 1 )V 1 C 20o -e o x(B pc C 20o -8 d C 0o )+^ BjC^ - x 2 B pc C 8 C 0 J (60) 

6$ equation: 

- £ A 2 [(H^) Opt »♦)']'- «♦" + - A l>{ ( 5 " 2 2 ~ ) [- S 29 0 + ( 9 pt + *> C 2e 0 3 

+v"w"[c 2?o + (e p t + 9) s 2 8 0 ] )- t - (kij - kSj) tSe 0 c 9o + <e pt + + + *)c 20o ] 

+2(km 2 S 0o v' + ^ 1 C 0o w’) +km 2 ^+^* [2 xC 0q (J + *) + xC^w’ - w] 

" ‘ W *Vi (61) 
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6* equation: 

Q#+ T ‘ C « S 0 O + * " 8 d c e 0 +*) ♦ *o jT 1 (vSe o+ * + wc 0o+ *)dx + i 5,0^*, 

+ 3 (6 pc + c. s e o+ *-e d C9 o ^)(;,c0 o ^ + e d S9 o ^) + i l6pc f l (vc 6o - ws 6o )dx 
+ (Bpc ♦ C«S0 o +* - 8 d Ce o+ 0) *(vC0 o+ * - CS 9o+< )dx 

♦Cc.Ce 0+ * + 6 d s e^*)^ 1 x(vs 6o+ 0 + iH: 9o+# )dx + s 2( e o+ * ) jf ' d* 

+ C 2 ( 0 O +*)^ w d * + (£ « 2 ' t «i>[ s e o C 0 o + *c 2eo + c 20o jT 1 ( 0 pt + *)dx] 

+2S PC So (v^ + 50 )dS + 26 j .£* X(c,v- 0 d 6 )d 5 

^(C^-BdCe^jr 1 (-4^) (*’♦• -f o'd' )ds + ^t« 2 S eo <r(l) + 2k* lCeo i(l) 

+2 r <* s e 0 + i *e 0 >/ 0 * (v’v' + 5 'w , )dx 1 dx^jJ 1 <vse 0 *«ce 0 >(c.$- fld 6)d* 

+ V ^ + 5.x)8 d*-^ 1 (C-8 d x)$ dJc + k* 2 Jf 1 0 + 

+2 * S(t, 9 -B d w)dx + » jT 1 (0 2 +w 2 )dx + k IB 2* 

*i(l6j c e 0 * + 2c e 0 £ 50 dx + c| o j T : xw' dx-^ 1 a dx>| JT 1 (v + ; .x)6 dx 

-s 0Ot £ l ^^5 l S)W 0o . % e o )dx-c l s eo (4ifi)^ i c 20o (5 + 5l) , dS 

- C 8o V (* 2 + 25 jx) (0 pt + ♦, (v + c 8 x)dx + S 0o C 0o f o ' X 2 (0 pt ♦0)(w- 8 d x)dx 
+C0 ° ^ ( v i 2 *x 2 — ) (w- 8 d x)dx-C 0o JT 1 x 2 (v + c,x)JT ?'w" dx! dx 

(jf 1 » dR-^).v lC2eo jT 1 x(0 pt + 0)(C- e d x)dx 

+, o ( 8pc c 20 o - 8 d c 0o ) (y+/‘ XV dx)-| e,c 0o (^ + jT 1 ** d x) 

~2 C loV TO ’(v+C,x)dx + C e o X(e 0 w‘ +vw’C 0o + 8,v)Cv+c 8 x)dx 
^OoBpcC. (t + ^ 1 525 d5 )*^ 1 x2 ^C 20o -wS 20o )dx-^ 

" 2 * 1 * (T + ^ 1 dx) + 0C 2 0 o JT 1 x 2 (e pt + 0)(O-6 d x)dx + 02c 2eo (jT 1 x 2 w dx-^) 

-**IKU 0 (J f 1 » dx-^)-2 ^ Ceo xv (5 - 6 d x)dx 

+(S e o ^ C 0 o > /a 1 ^eptMW-BdXjdS-ViC^^ 1 (e pt+ 0 ) 7 (w-B d x)dx 

+9 i C 0 o fo 1 (v+4 8 x) 5 dx- (S 0o+0 C 6o ) JT» x(v + c,x )0 dx 

~ 2 C *ofo' *^ 0 Pt ♦ 0 )^(v + C 8 x)dx - 2 (S 0o+0 C 0o ) X 1 x6(«-e d x)dx 

" C 0 O fo 5 < e pt + *) 5 (O-e d x)dx + 2v i C 0o jT 1 5 (w - 8 d x)dx 

+C0 o ^ (*♦»!>(»♦«,«>* dx-f C 0o / o J x(v + C 8 x)0 dx-f c e 0 *{/ 0 1 5v dx + ^) 
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Terms not satisfying equation (3) are 

-2(^ + i)(km 2 SQ o 6v + km l C 0o 6w)|J 


- v i 

Bo “ U “ <R)6o 

(63) 

(64) 


The terms of equation (64) must be included in the analysis as discussed above 
in order to obtain the correct results. Equations (59)- (62) are the hybrid 
nonlinear equations of motion that will be solved b> Galerkin's method. 


Application of Galerkin's Method 


In transforming equations (59)- (62) into modal equations we use the mode 
shapes for a nonrotating cantilever beam. We also assume that the motion is 
characterized by small perturbation motions about a steady equilibrium opera- 
ting condition that depend on dimensionless time i|/(-8t). Thus, 


N 

v - E [V Q i + AV 1 (t|/)]f 1 (x) 
i-1 

N 

w - 2 t w oi + AW i (ij»)]'l , i (x) 
i-1 

N 

♦ ■ J2 [*oi + ] 0 1 (x) 


(65) 


where 


4» i (x) - cosh(B^x) - 008(8^^5) + a^[sinh(6ix) - sin(B i x) ] 
0 1 (x) - & sin(y i x) 


(66) 


The constants aj_ and 8^ are defined in reference 6 and y^ - rr [ i — (1/2)]. 
We also assume that 


and that 


* - *o 



+ 

(67) 

- 0 t* 

(68) 
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Nonlinear algebraic equations for V D j , W 0 j, t 0 j , and determine the 

equilibrium deflections. The perturbation deflections are governed by a set 
of homogeneous ordinary differential equations for AVj, AWj, A$j , and A$, 

with constant coefficients depending on V 0 j, W Q j , 4> 0 j , and * c . The modal 
equilibrium equations are as follows: 

<5V ol equation: 

2 ^A 2 Bj 4 « i jV 0 j - (A 2 - Ai) j^S 0o 8j 4 6 1 j - S 2 q o *ok v klj " e t p ij) 

^ 20 o E ♦ok [ E WkUj-ZetViijjJv^ + fAz-A!)^ -y* Bj^j 

+C 20 O ( E *ok v kij " 9 t p ij) + S 20 o E *ok ( E ♦o£ Y k£ij ” 2e t v kij)l w oj 
'k-1 * k»l >£■! /J 

+(Mij +5 1 L i j)V 0 j " V o j6ijC c 0 o - *q S 20 o - ♦o 2c 20 o ) +W oj fi ij( S 0 o C 0 o + *o C 20 o - ♦o 2s 20 o > 
+Aj ji 1 B p c* o C0 o + e o * o ^S0 o + y 2 ^ + B j*o (0 2 + ♦ 0 5 2 > 6 1J 

+ e { ^iC 26o Q i: j - Se^e^j)^ - * o [ s 20 o c j + c 20 o (*o c j “ e t D j " 2 ^i B j>J 6 ij"*o R ij*oj|) 

" l*o C Q 0 - ( B pc S 0 o + C s> l A i • B 1 S 0 O B 1 + 6 [( c i + 2 ® l B l) s 0 o ” 0 t c i s 6 o c e o 
■^i^oAi-Vi^eoBl-Q^e^!)-^ <*] 1-1,2 N (69) 
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6W oi equation: 


£<(A 2 -A,)[- ^2 Wkij * e t p ij) 

+S 25 0 Z» *ok^Z) *ot T klij ■ 2e t v ilJ^j v o] +A lSj'*«ijW 0 j+ (A 2 - Al)f8; o fj 4 <ij 

“ s 2§ 0 ( 2 Z *ok v kij ” 0 t p ljW 2 g £ *ok ( E *o* Y kUj " 20 t v klj) w oJ 
\k-l / k-1 \&-l U 

+ (MiJ ♦ *iWj)w 0 j + VojSij (Se o c 0o + ♦ o c 20o - * o 2 s 20o ) - (s^ o + V2e o + *o 2 c 2 e 0 > 


+A J [“ 5 l 0 pc*o S e o + 5 oV c 0 o “ *O S 0 O > l 6 lj + B J*o<*2 ■ *o B 2> 6 lj 

+ s{- (R iJ +25 l c >lj )c e 0 * 0 j + ]C u ljk v oj w ok c e 0 + (*o - f c «o) c 8 o 0 1 J w oj 

+ 8l c 9o I « V oj - ‘o C 2e 0 <Cj + 2S I » J )« ij J 

t«l (B pc C 6o - B d ) + 5 o S 0o ]A 1 - BiCe^t +* [(C t + 2S 1 B 1 )3 9 o C 9o - e t C 9o (D 1 + 2iiC t ) 

"^l C 20 c (®i + ®l A l) " ®O^ 0 pc C 20 o ” B d C 0 o ^ B l + 0 1 \ C 0 o B i " ^pC C 8 C 0 o C lJ 
64> 0l equation: 

21 + KY j 2 6i j H 0 j + ^ A 2 “ A i ) | J) v ijk£" J~ ^ W oj W ok “ V oj V ok^ 

] N 

+ 0 t r v; jk [c 20o (w oj w ok -v oj v ok )+s 20o v oj w ok ] 

N N 1 

+ E E ’fijkl*o J [ C 29 0 «okWot-Woll) + S25 0 Vok w ol] 

1-1 k-l f 

+(iS 2 +e j *„)« 1j c 29o+ ^ c| 0 Q^w OJ ) 

■ ^ < Vi F i ' ( kJj ,- Ei I )< s 9 o c 9 o Ei - e t F iCa 9 o )+ V 9 t F i <»> 
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equation: 
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The linearized perturbation equations are easily expressed In matrix form as 


[M]{AX} + [C] {AX} + [K] {AX} - 0 


(73) 


where [M] , [C], and [K] are, respectively, the mass, gyroscopic and damping, 
and stiffness matrices given below and 


[M] 



{AX} - - 

A Vl j 

AWi , 

|A#i 


(74) 

{ U ! 

0 

\ 

o I 

1 

-(W or 6 d Bi) 


t 

0 1 

S ij( 1 + 2f) ' ' 

1 

_ ill s 

96 S Ji 1 

^ s ij 1 
— — — •— ■! J 

^n 2<S ij | 

Vol+CeBi-^ At 
^m 2 ®i 

(75) 

-(W ol -B<,Bi) 
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1 

L 1 
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N> 

00 




where the node shape Integrals are 

*i - C 'i dS 

B t - Jf > «* t dS 

C 1 ‘ f a ' x 2 * 1 dx 

D 1 " /q x 3 * i dX 

E 1 “ e i d * 

F i " f Q l dx 

G i * f 0 l x 2 ©i dx 
Hi - f Q l x 3 0i dx 

l a - V »i*j d * 

J lj - f Q l d * 

K ij ' V 50 i 0 j d * 

L U - to < 1 -S>W d * 


«1J vr dS 



°ij - **iV d * 


P U * d * 

q ij * 4 1 5y t°j d x 

oij • * 0 iY d * 

R 1J ■ f 0 l * 2 1'iSj dx 
S H - •£' ’ r l 0 j d * 
s ij ‘ f a ' 0 i y .i ' d * 

T ljk - -^TT fo VVV diE 

u ijk - f a ' « 

-jf 1 * 2 »i C vv d *i di 

V ijk - f 0 ' ¥)"V dx 
v ljk - f a ' 50 i T j" y k" dx 
“ijk - / 0 ‘ Iffy dx 
X ljk -C ™ffy d * 

Y ijkl- f 0 ' 'Wk'V’ dx 

Z iJ * f 0 ' *'fi ’ d * 

-f 0 l * 3 f* dx, dx 


( 77 ) 
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Many of these integrals have been evaluated in closed form by use of 
references 7 and 8. Additionally, these and all the remaining integrals were 
evaluated numerically. The matrix [K] is simply the Jacobian of equa- 
tions (69)- (72) and will not be written out in detail. Representing 
equations (69)- (72) as 


Then 


where 


Y t - 0 i - 1 , 2 , . . . , 3N + 1 


(78) 



(79) 


<V)> - l'' 0 j.w 0 j ,* 0 j .* 0 J T 

We note that [M] and [K] are symmetric and [C] is antisymmetric for 
Y - 0 (in vacuo). Thus, equation (73), a standard eigenvalue problem, governs 
the stability of small motions about the equilibrium operating condition. 

As described above, it is possible to use the torsion moment boundary 
condition, equation (6), instead of the root torsion equation (62). If this 
is done, the dimensionless form of equation (6) becomes 

k A 2 . 

Q* - *c*'(0)+-y- [♦'(0) + e pt (0)] 

+Z b {A 1 w"(0) - (A 2 - A!)S 3 o [*"(O)C 0 o - w M (0)SS o ]} 

+8 d {A 2 v"(0) - (A 2 - A 1 )Se o [V"(O)Sg o + w"(O)C 0 o ]} (80) 

The modal equilibrium equation is 



N 

- <*2 - M)Se (Vote, - w o1 s, ) ] 

i-1 

N 

+28 d 2 6i 2 [A 2 V oi - (A 2 -Ai)Se (V ol S^ o + W oi Cg o )] (81) 

i-1 


The perturbation equations may be easily solved for and substituted into 

the other equations since there are no time derivations. Thus, 
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A* 






£ st^ w i 

i-l 


-*-2B d A 2 2 •i^-ZC.Ma-AOSS £ B 1 2 (AV 1 Ce - AWiS§ ) 
i-l ° i-l ° 

N 1 

-2B d (A 2 -A 1 )Se o 2 3 i 2 ^ v i s 0 o + AW i c 5 o )J 


( 82 ) 


Numerical results obtained with equations (81) and (82) are virtually identi- 
cal with those using the 5$> 0 and 6A$ equations derived above based on inte- 
grated torsion moments. The use of equations (81) and (82) is considerably 
simpler and provides a reasonable check for numerical results. 


Modal Analysis 

We now describe a modal analysis that greatly simplifies numerical 
computation. From equation (73) the stability of small motions about the 
equilibrium operating condition is determined by the eigenvalues of the 
6N + 2 x 6N + 2 matrix [P] where 

li-[ ° 1 1| AX }- [ p i | AX ) < 83 > 

l ax) L-m-'k -M“*C J Jax / (ax) 

Since we are primarily concerned with lower frequency instabilities 
(first lead-lag, first flap, and first torsion frequencies), there is a value 
of N for which any increase in N will not appreciably change the eigen- 
values associated with these lower frequencies. It is at this value of N 
that the eigenvalues are considered to be converged. For practical hingeless 
rotor configurations, N - 5 gives suitably converged eigenvalues; the matrix 
[P] is thus 32 x 32. By a change of modal coordinates, the size of the matrix 
[P] may be greatly reduced without significantly changing the eigenvalues of 
interest. Such a transformation may be found by first considering free vibra- 
tions in vacuo of the blade about the equilibrium deflected state. The equa- 
tion of motion, analogous to equation (73), is 

[M S ](AX} + [G]{AX> + [K V ]{AX' - 0 (84) 

where the subscripts s and v imply the symmetric part and the vacuum case, 
respectively. Both [M s ] and [Ky] are thus symmetric; [Ky] is equal to [K] 
with all the aerodynamic terms set equal to zero. The matrix [G] is anti- 
symmetric and equal to [C] with all aerodynamic terms set equal to zero. The 
presence of [G] causes the eigenvectors of free vibration to be complex. This 
may be avoided for computational efficiency by approximating equation (84) as 

[M s ] {AX} + [Ky]{AX> - 0 (85) 
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The matrix of the eigenvectors [U] is orthogonal with respect to [Mg]. It is, 
therefore, approximately orthogonal with respect to [M] as well since the 
elements of the antisymmetric part of [M] are very small. Thus, 

[U T M,U] - [I] * [U t MU] (86) 

According to Melrovltch (ref. 9) a so-called principal coordinate transforma- 
tion for equation (73) may be determined by replacing (AX} by [U](AX>. We may 
then premultiply equation (73) by [Ul T to take advantage of the form of 
equation (86) yielding 

[I] {AX} + [U T CU](AX) + [U T KU]{AX} - 0 (87) 

* 

Hence, 



The matrices [P] and [P*] have virtually the same eigenvalues. However, 
because of the nature of this modal coordinate transformation from [P] to [P*], 
the rows and columns corresponding to nigh frequency modes of both [U*KU] and 
[U T CU] may be removed without affecting the eigenvalues of the low frequency 
modes of Interest. These 3N + 2 * 3N + 2 matrices are thus reduced to M * M 
matrices whore rows and columns correspond to the M low frequency modes that 
are retained. The rows and columns that are retained In [U T KU] and [UTCU] may- 
be chosen in two ways: (1) the M rows and columns that correspond to the M 
lowest frequency modes of the blade may be retained, or (2) the M rows and 
columns that correspond to M modes selected arbitrarily from the lowest 
lead-lag, the lowest flap, and the lowest torsion frequency modes are retained. 
For the second case, under certain conditions, M • 4 or 5 will result in con- 
verged eigenvalues. In either case, suitably converged results do not require 
M > 8. 

The reduced matrices are analogous to stiffness and damping matrices 
generated from M coupled, rotating modes. Since the analysis is formulated 
in terms of standard cantilever mode shapes, however, repeated numerical inte- 
gration of modal Integrals is not necessary for different values of blade 
stiffnesses. Instead, the matrix operations described above lead to a net 
savings in CPU time. 


CONCLUDING REMARKS 


Hybrid equations of motion are developed for an elastic blade canti- 
levered in bending and having a torsional root spring to simulate pitch link 
flexibility. The blade is assumed to have coincident mass center, tension 
center, aerodynamic center, and elastic axes. Droop, precone, twist, sweep. 
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torque offset, and blade root offset are included in the model. Quasi-steady 
aerodynamic loading la assumed to be adequate to investigate the low frequency 
type of unstable motion common in hingeless rotor systems. The solution is 
obtained by Galerkln's method and a modal analysis. The stability of small 
motloi.s about the equilibrium operating condition is governed hy a standard 
eigenvalue problem where the elements of the stability matrix depend on the 
solution of the equilibrium equations. In the analysis, two different forms 
of the root torsion equation are developed. One is based on the torsion 
moment boundary condition at the root of the blade and the other is based on 
integrated torsion moments derived from the kinetic energy. Numerical results 
for the two cases are virtually identical providing a reasonable check of the 
equations. 
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